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Concurrent remote entanglement of distant, non-interacting quantum entities is a crucial function 
for quantum information processing. In contrast with the existing protocols which employ addition 
of signals to generate entanglement between two remote qubits, the continuous variable protocol we 
present is based on multiplieation of signals. This protocol can be straightforwardly implemented 
by a novel Josephson junction mixing circuit. Our scheme would be able to generate provable 
entanglement even in presence of practical imperfections: finite quantum efficiency of detectors and 
undesired photon loss in current state-of-the-art devices. 


Generation of entangled states between spatially sep¬ 
arated non-interacting quantum systems is an indispens¬ 
able ingredient for large-scale quantum information pro¬ 
cessing W- In particular, concurrent remote entan¬ 
glement, in which propagating quantum signals do not 
interact with both the systems under consideration, is a 
desirable feature of a scalable module-based architecture 

[SHE]. 

It is well known that a non-linear operation is nec¬ 
essary to achieve entanglement. Existing protocols for 
heralded concurrent remote entanglement employ linear 
optical elements in the processing stage (e.g. beam split¬ 
ters), while the necessary nonlinearity is provided in the 
final stage by photon detection [9HI5]. While these meth¬ 
ods relied on addition of signals using beam-splitters to 
erase ‘which qubit’ information, our proposed method re¬ 
lies on multiplication of signals coming from each qubit 
to delete their local orientation. This multiplication is 
achieved by a new type of nonlinear signal processing. 
Josephson junction based superconducting circuit QED 
systems have access to strong, tunable, purely dispersive 
nonlinearities, making them natural candidates for im¬ 
plementing this protocol. In fact, sequential remote en¬ 
tanglement with linear microwave signal processing has 
already been performed using Josephson junction circuits 

m- 

The first step of our protocol generates local entan¬ 
glement na HH] between a stationary superconducting 
qubit (for definiteness, a transmon qubit) and a prop¬ 
agating microwave mode (cf. (I) in Eig. for each 
of Alice and Bob [19]. Both Alice and Bob are ini¬ 
tialized, using local 7r/2 rotations (Y^/^), to a super¬ 
position of their ground (|^)) and excited (|e)) states, 
given by: (|^) + |e))/\/2. Propagating modes, with 
coherent states of amplitude and /3o, and tempo¬ 
ral profile and cos(cjf^t)0(—t), 

are incident resonantly on two cavities, exciting their 
fundamental modes A and B, with frequencies (decay 
rates) uJa(f^A) and with Ka(b) ^ ^a{b)‘ These 

modes interact dispersively through cross-Kerr interac¬ 
tion Eniin] with Alice and Bob. This operation is re¬ 
ferred to as the conditional displacement gate (CD). It 
imparts a qubit-state-dependent phase-shift on the out- 





FIG. 1. (color online) Remote entanglement protocol 
schematic. The first step of the protocol (I) consists of en¬ 
tangling two stationary qubits, Alice (in dark red) and Bob 
(in dark green), with two propagating modes (indicated by 
thicker lines), ancilla signal a and ancilla signal b respectively, 
each initially prepared in a coherent state. This is achieved 
by first applying a 7r/2 rotation (Y^^^) on Alice (Bob), fol¬ 
lowed by a conditional displacement gate (CD) on Alice (Bob) 
and ancilla signal a (ancilla signal b). In the next step (II), 
a nonlinear interaction of the ancilla signal a, ancilla signal 
b and parity signal c (implemented by the Josephson Para¬ 
metric Multiplier (JPM), see below), followed by a homodyne 
detection of the parity signal mode, effectively realizes a joint 
two-qubit parity measurement. The resulting qubit-photon 
state has either even or odd joint qubit parity conditioned 
on the integrated homodyne current being in either lobe of 
the distribution (indicated by dashed lines). The last step 
(III), which disentangles the qubits from the photon states, 
comprises of homodyne measurements of the a and b modes, 
denoted by HDa and HDb- Conditioned on the measurement 
outcome in (II) and (HI), the qubits are projected onto the 
even or odd Bell manifold, with a relative phase that depends 
on the measurement outcome in (III). 

going microwave modes. The resultant entangled qubit- 
photon states output from Alice and Bob’s cavities can be 
written as: (|e, a)p\g, —a))/\/2 and (|e, /d) + |^, —P))/\/2 
[22] with temporal profiles cos(cc;at)0(—t) and 

cos(co’5t)0(_t), respectively [SSlISJ. Without loss 
of generality, we may assume a, /3 G and they need not 
be equal in our protocol. 

In the next step, we realize a joint two-qubit parity 
measurement by first capturing the propagating modes 
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in resonators and then employing a nonlinear dissipation 
process. To that end, we introduce the Josephson Para¬ 
metric Multiplier (JPM) (see step (II) in Fig. [^. The 
JPM comprises three resonators and a nonlinear four 
wave mixing element, the Josephson Four Wave Mixer 
(JFWM) (Fig. [^. The three resonators have fundamen¬ 
tal modes (frequencies, decay rates) 3.{uJa, i^b) 

and c{ujc^ i^c)- The outputs of cavity modes A and B, af¬ 
ter propagating through transmission lines, act as inputs 
to the a and b modes, respectively. Due to their partic¬ 
ular temporal profiles, these flying modes are perfectly 
captured at t = 0. A coupled two-mode dissipation is 
then turned on at t = 0, which removes pairs of photons 
from the a and b modes at a rate /^ 2 ph- This dissipa¬ 
tion, mediated by the jump operator ab, is realized by 
the JFWM, together with the dissipation of the c mode 
in the following way. 

The JFWM consists of four nominally identical Joseph¬ 
son junctions, as shown in Fig. [^and has four interact¬ 
ing normal modes, which are negligibly shifted in fre¬ 
quency from the original modes a, b, c, in the presence 
of a stiff, off-resonant pump mode with frequency cho¬ 
sen to be uOp = ujc — — ^b- Thus, under the rotating 

wave approximation, the mode-mixing arising out of the 
Josephson nonlinearity leads to an interaction Hamilto¬ 
nian of the form = ige~^^p^abc^ -f- h.c., where g, 

the effective interaction strength, depends on the pump 
amplitude [24]. If the cavities are designed and pump 
strength is chosen such that 

f^a.f^b < g,l^2ph < ( 1 ) 

the JPM will provide unidirectional conversion: photons 
in modes a and b are converted into the c mode, which 
leaks out before it can be converted back into the a and 
b, providing the desired two-photon dissipation channel, 
characterized by the decay rate /^ 2 ph = ^9^ j after adia¬ 
batic elimination [25]. The nonlinear dissipation channel 
is monitored with a homodyne detection scheme, denoted 
by HDc, with phase angle arg(g(a/d), which measures the 
value of the integrated homodyne current Xc- By select¬ 
ing outcomes Xc, which will fall on either lobe of the dis¬ 
tribution, centered at ±.2gmhi{a^ shown in (II) in 

Fig-lU the qubit-photon state is projected on to the even 
or odd joint qubit-parity subspace. The error rate at this 
step can be controlled by selecting extremal outcomes, 
beyond some cut-off in each lobe of the distribution. At 
the end of the measurement, the two-photon dissipation 
is turned off by switching off the pump at ujp. 

The homodyne measurement of c in step (II), is fol¬ 
lowed in step (III) (Fig.[^ by homodyne measurements of 
the modes a and b, denoted by HDa and HDb- This last 
pair of measurements is crucial because, while the two¬ 
mode dissipation projects onto the even or odd qubit- 
parity subspace, the photons left over in the modes a 
and b after step (II) are in a two-mode squeezed state 
which remains entangled with the state of the qubits. 
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FIG. 2. (color online) Schematic of the JPM of Fig. 

The modes a(ct;a), b(ct;6) and c{uc) when pumped at Up — 
uJc—^a—^b (orange) participate in a non-linear, three-wave in¬ 
teraction Hint/^ = -h h.c.. This nonlinear mode 

mixing arises out of the Josephson Four Wave Mixer (JFWM). 
(Inset) The JFWM has four nominally identical Josephson 
junctions connected electrically as shown. The coupled sys¬ 
tem has four mutually orthogonal normal (electrical) modes 
[24] , which couple non-linearly, corresponding to the cavity 
modes a, b, c and the pump mode. 


Step (III) disentangles the qubits from these microwave 
modes, as follows. 

Consider the case in which the two-qubit parity mea¬ 
surement projects the system to the even two-qubit par¬ 
ity subspace. The sign of the X quadrature measure¬ 
ments Xa^xi) is correlated with the probability that the 
qubits are in the \gg){xa^xi) < 0) or |ee)(xa,X 5 > 0) 
states and only in certain regions of the (xa^x^) plane 
along the line Xa = —xi) are the two qubit states strongly 
entangled (Fig[^ upper panels). Conversely, Y quadra¬ 
ture measurements will give results centered around ya = 
^5 = 0, are not correlated with the two qubit states, and 
do not distinguish between them; the result is that Y 
measurements always entangle the two qubits (with a 
relative phase which interpolates between the even and 
odd Bell states) (Fig. lower panels). Similar reasoning 
holds for the odd two-qubit parity subspace outcomes of 
step (II). Hence, while it is possible to have reasonable 
success rate by making X measurements on modes a and 
b, it is always preferable to measure the Y quadrature 
for optimal success rate. 

While it is possible to perform complete stochastic 
master equation simulations of the protocol we have just 
outlined [24] , given the assumed separation of time scales 
(Eqn. 0 ) , we have used in Fig. |^a simpler approximate, 
but accurate, model, which can be solved analytically and 
provides physical insight. Following the capture of the 
propagating microwave modes in signal resonators, the 
state of the system comprised of Alice, Bob and modes 
a, b, at t = 0, is given by p(t = 0) = |'0)('0|, where = 
{\ee,a,j3) + \gg,-a,-j3) + \eg,a,-l3) + \ge,-a,p))/2. 
Since the two-qubit parity measurement of step (H) de- 
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FIG. 3. (color online) Probability distribution of outcomes, 
overlap (F) with Bell-state = (|ee) + \gg))/V2, concur¬ 
rence (C) of the joint two-qubit system pq and gradient of 
overlap are plotted when the qubit-photon state is projected 
onto one with even two-qubit parity after the second step of 
the protocol. The top (bottom) row panels correspond to ho¬ 
modyne detection of the X (Y) quadratures of modes a, b. 
We choose a = /3 = 0.75 and assume perfect quantum effi¬ 
ciency and zero spurious photon loss. Panel (a) shows the 
probability of outcomes P{xa,Xb) for X measurements, corre¬ 
sponding to which, we see that the two-qubit state is projected 
on to 1$^) for values around the diagonal Xb = —Xa (panel 
(b)). For events occurring in the quadrant Xa,Xb > 0(< 0), 
the two-qubit state is projected on to \ee)(\gg)). Panel (c) 
shows the concurrence C{xa,Xb) of pq for the different out¬ 
comes, which varies from 0 (for pq being a separable state 
|ee) or \gg)) to 1 (in case of maximum entanglement). Panel 
(d) shows the gradient of the overlap as a function of Xa,Xb, 
which is zero (in green) for a narrow region around the line 
Xa = —Xb where an entangled state is obtained. It changes 
rapidly on moving away from the line Xa = —Xb and goes 
back to zero when the qubit state is projected on to |ee) or 
\gg). Panel (e) shows the probability of outcomes P(ya,yb) 
for Y measurements. For events around the line ya = —2/6, 
the two-qubit state is once again projected on to |^^). How¬ 
ever, for Y measurements, the phase of the generated Bell 
state varies continuously, depending on the particular out¬ 
come (//a, 2/6), indicated by the existence of alternating bright 
and dark fringes in fidelity (panel (f)). For instance, along 
the line 2/a = 2/fe? fhe two-qubit state oscillates continuously 
between and |^~). Panel (g) shows the concurrence 

C{ya,yb) of Pq, which is 1 for all measurement outcomes, in¬ 
dicating generation of maximal entanglement for all outcomes 
(2/a,2/fe)- The small regions of low entanglement are artifacts 
of the analytic approximation which replaces stochastic evo¬ 
lution of the system with a deterministic Lindblad evolution 
(see further discussion in the text). Panel (h) shows the gra¬ 
dient of overlap which is zero where the two-qubit state is 
projected on to and changes rapidly as the two-qubit 

state oscillates between and |^~). 





pends only on unambiguously inferring which side of the 
distribution the outcome is on, for the purposes of ana¬ 
lytic computation, we separately average over the differ¬ 


ent outcomes for the two lobes of the distribution (see 
Fig. [^. This amounts to replacing the stochastic evolu¬ 
tion of the whole system by separate deterministic Lind¬ 
blad evolutions for the even and odd qubit parity sub¬ 
spaces. During this evolution, the single photon losses of 
modes a, b are negligible due to Eqn. 0. The system 
density-matrix in even (odd) qubit-parity subspace Pe(o) 
thus evolves according to: 


= K2phr>(ab)/)e(<,), (2) 

where /5e(o)(i = 0) = \'ipe(o)){i’eio)\,\'>Pe) = {\ee,a,(3) + 
\gg,-a,-/3))/V2, lipo) = {\eg,a,-(3) + \ge,-a,(3))/V2 

and P(0)pe(o) = Ope{o)0^ - (0’*‘Ope(o) + Pe(o)0’*‘0)/2 

is the Lindblad dissipation operator. The quasi-steady 
state at the end of this evolution, denoted by [24] . 
subsequently evolves under the single photon loss of a and 
b that are monitored by HDa and HDb- The resulting 
measurement of both S G {X, Y } quadratures of a and b 
modes results in the system density matrix evolving to: 


qs 

Pe{o) 






Ms 


|Ca,6)(Ca,6|- (3) 


The post-measurement two-qubit density matrix pq is 
computed by tracing out the modes a and b. Conditioned 
on the outcomes in (II) and (III), the qubits are pro¬ 
jected onto an entangled state in the subspace spanned 
by {\gg)^ {\^9)i The continuous nature of 

entanglement generation appears as a relative complex 
amplitude of the two terms of the Bell state, which is de¬ 
termined by the measurement outcome in HDa and HDb- 
In Fig. we show the probability of outcomes, overlap 
with the Bell-state |4>T) = (|ee) + concurrence 

and gradient of the overlap for either X measurements 
(top row panels) or Y measurements (bottom row panels) 
of the a and b modes. For X measurements, we see that 
the majority of the events occur for either Xa^x^ > 0 or 
Xa^Xb < 0, hence projecting the qubit state onto product 
states |ee) or \gg). However, the (non-negligible) number 
of outcomes near the line Xa = —x^^ do project the qubit 
onto the entangle state Accordingly, the concur¬ 

rence C{xa^Xb) varies from 0 (for pq separable) to 1 (in 
case of maximum entanglement). The width of the region 
in phase-space where entanglement is generated is a func¬ 
tion of a, /3 and decreases as a, /3 are increased. The rate 
of variation of entanglement is indicated by gradient of 
the overlap |<I>T) which varies most rapidly perpendicular 
to the line Xa = —Xb- For Y measurements,in contrast,a 
maximally entangled state is generated for all outcomes, 
with the phase of the generated Bell state varying contin¬ 
uously in the form = (|ee)+e*'^|^^))/v^, giving con¬ 
currence equal to unity at all points. An increase in a, /3 
makes this variation more rapid. In this case the gradi¬ 
ent of overlap is not a measure of entanglement, but just 
describes the variation of the phase, ip{ya,yb)- A similar 
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FIG. 4. (color online) Maximum overlap with Bell state 
is shown upon variation of incident state amplitude a 
and efficiency parameter 77 for measurement of Y quadratures 
of a and b. A sample of 500 trajectories were simulated for 
each data point in the (a, 77 ) space and maximum fidelity was 
noted. For perfect efficiency (77 = 1), it is always possible to 
generate |4>+). The maximum target fidelity goes down as 77 
is lowered. Higher values of a is more susceptible to finite 
efficiency. However, even for 77 = 0.7 and a = 0.5, we obtain 
fidelities in excess of 80%, indicating the robustness of the 
scheme to photon loss and finite quantum efficiency. 


computation for the odd manifold shows similar results 
for X and Y measurements, with |ee) ^ \eg), \gg) \ge) 
with the features in the Fig. [^rotated by 7 r/ 2 . 

As mentioned before, our simplified analytic model 
does not track the precise value of the homodyne cur¬ 
rent Xc, but only its sign. This corresponds to averag¬ 
ing over different outcomes with even (odd) qubit parity, 
leading to Eqn. <© making the quasi-stationary state 
slightly impure. This leads to impurity in the post¬ 
measurement qubit state for a some regions in the out¬ 
come plane and is therefore an artifact of the approxi¬ 
mation. The small regions of spurious zero concurrence 
in Fig. are due to this. As a, /3 are increased, the 
impurity due to the Lindblad evolution increases. This 
restricts the accuracy of the approximate analytical the¬ 
ory to P ^ 1. For P ^ 1, one can numerically 
simulate the evolution using the stochastic master equa¬ 
tion. A comparison between the analytical solutions and 
stochastic master equation solutions is provided in [24] 
for 13 

In what follows, we test the robustness of our protocol 
to imperfections arising out of undesired photon loss and 
finite quantum efficiency. We treat these imperfections 
together as a general efficiency parameter 77 . We present 
here results of numerical simulations for the case when 
both Y quadratures were measured (Fig. |^. Similar re¬ 
sults are obtained when X measurements are performed. 


While for perfect efficiency (77 = 1), it is always possible 
to generate an entangled Bell state, even for 77 = 0.7 and 
a = 0.5, fidelity in excess of 80% is obtained, indicat¬ 
ing the robustness of the protocol to these imperfections. 
However, the photon losses before the JPM prevent a 
complete trade-off between the success-rate and fidelity 
as will be discussed in [26]. Note that a low value of a, /3 
lowers the success rate of entanglement generation since 
the two-qubit parity measurement in the step (II) relies 
on unambiguously inferring the location of the homodyne 
outcome of the parity signal. 

To summarize, we have presented a protocol for re¬ 
motely entangling two qubits by performing a set of con¬ 
current quantum operations on propagating microwave 
modes entangled with the qubits. In contrast to exist¬ 
ing schemes based on linear optical elements and photon 
detectors, we propose a qualitatively different approach, 
based on the multiplication of quantum signals prior to 
continuous measurements to generate remote entangle¬ 
ment. This multiplication is achieved using the Joseph- 
son nonlinearity and the high detection efficiency of mi¬ 
crowave radiation in circuit QED systems promise a much 
higher success rate of entanglement generation compared 
to its optical counterparts. 
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